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Abstract: We compute the electroweak (EW) contributions to squark-squark pair pro- 
duction processes at the LHC within the framework of the Minimal Supersymmetric Stan- 
dard Model (MSSM). Both tree-level EW contributions, of ©(a^a -|- a^), and next-to- 
leading order (NLO) EW corrections, of 0{a1a), are calculated. Depending on the fla- 
vor and chirality of the produced quarks, many interferences between EW-mediated and 
QCD-mediated diagrams give non-zero contributions at tree-level and NLO. We discuss 
the computational techniques and present an extensive numerical analysis for inclusive 
squark-squark production as well as for subsets and single processes. While the tree-level 
EW contributions to the integrated cross sections can reach the 20% level, the NLO EW 
corrections typically lower the LO prediction by a few percent. 



Keywords: 



Bupersymmetry Phenomenology, NLO Computations, Hadronic Colliders . 



Contents 



|l]. Introduction 



|2|. Classification of processes and tree-level cross sections 



i 



i 



Virtual and real corrections of 0{a.^a) 

Virtual corrections 
Real corrections 

Real photon emission 



3.2 



3.2.1 



3.2.2| Real gluon emission 



3.3 



3.2.3 Real quark emission 

Factorization of initial-state collinear singularities 



Numerical results 
[4.1| Input Parameters 

Total hadronic cross sections 
Differential distributions 



4.2 



El 



[sj. Conclusions 

|A] . Tree-level cross sections 



Bremsstrahlung integrals 

Soft and collinear photon bremsstrahlung 



B.l 



B.2 Soft and collinear gluon bremsstrahlung 
|B.3| Collinear quark bremsstrahlung 







10 

|n 

13 
15 

21 

23 

26 

27 

m 



in. Feynman diagrams 



29 



1. Introduction 

Supersymmetry (SUSY) [1] is one of the most appealing scenarios for physics beyond the 
Standard Model (SM). Being proposed as the only nontrivial extension to the space-time 
symmetry, it was found to have many nice phenomenological properties which overcome 
weaknesses of the SM. With SUSY particles at the TeV scale or below, the technical part 
of the hierarchy problem is solved by stabilizing the electroweak scale, in particular the 
mass of the Higgs boson. In addition, many super symmetric extensions of the SM offer a 
dark matter candidate that gives the observed dark matter relict density [2] . Of particular 
interest is the Minimal Super symmetric Standard Model (MSSM) [3-5] which provides a 
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fit to electrowcak precision observables and B-physics data with a comparable to the 
SM and a naturally light Higgs boson [6-8] and can explain the measured value of the 
anomalous magnetic moment [9,10]. 

If supersymmetry is realized at the TeV scale, it will be probed at the Large Hadron 
Collider (LHC) at CERN. Indeed, the 95% C.L. area of the (mo,mi/2) plane of the Con- 
strained MSSM (CMSSM), lies largely within the region that can be explored within Ifb""*" 
of integrated luminosity [11]. Within the framework of the MSSM, imposing R-parity 
conservation, SUSY particles can only be produced in pairs. Among the potential SUSY 
discovery channels, the direct production of pairs of color-charged SUSY particles is of par- 
ticular importance at hadron colliders since it proceeds via the strong interaction. Many 
searches for squarks and gluinos have thus already been performed at high-energy colliders. 
Results of the D0 and CDF collaboration can be found in e.g. [12]. Studies for the LHC 
are based on Monte Carlo simulations. It has been shown that there is the possibility of 
early SUSY discovery within lfb~^ of data in the inclusive jets plus missing energy channel, 
provided that the SUSY particles are not too heavy [13]. 

First theoretical cross section predictions for squark and gluino pair production pro- 
cesses based on leading order (LO) calculations were made already many years ago [14-17]. 
Later calculations of next-to-leading order (NLO) in perturbative QCD [18, 19] could re- 
duce theoretical uncertainties considerably and revealed corrections of typically 20 — 30%. 
Recently also results beyond the one-loop level in QCD have become available [20-23], in- 
creasing the cross section by another 2 — 10%, and stabilizing the prediction considerably. 

For a reliable cross section prediction, also electrowcak (EW) contributions have to 
be taken into account which are of the same size as the NNLO QCD contributions. The 
contributing processes are manifold and their interplay is nontrivial, in particular if not 
only QCD-mediated but also EW-mediated production channels exist at tree-level. 

The latter arise from qq annihilation or qq scattering, as well as photon-induced pro- 
cesses, and contribute at 0{asa + a'^). For squark (anti)-squark production processes, the 
qq (gg)-channels can rise the LO cross section by up to 20% [24,25]. The photon-induced 
channels are typically more important for the pair production of lighter stops [26], and 
have a reduced impact on squark-anti-squark [27] and gluino-squark production [28]. 

NLO EW corrections contribute at 0{a^a) and have been investigated for stop-anti- 
stop [26,29], squark-anti-squark [27], gluino-squark [28] and gluino-gluino [30] production 
processes. In this paper, we provide the yet-missing NLO EW corrections to squark-squark 
production. We give details on the NLO computation and present an elaborate numerical 
analysis of all squark-squark production processes, including the anti-particles. 

In the context of all squark and gluino production processes, squark-squark production 
is of particular interest at the proton-proton collider LHC. The partonic process proceeds 
at LO from gg-induced diagrams only. Squark-anti-squark and gluino-gluino production 
require qq or gg initial states instead. Since the final-state SUSY particles are very massive, 
an important contribution to the hadronic cross sections arises from the high-x region where 
valence-quark densities dominate. As a result, squark-squark production has generally 
a higher tree-level yield than squark-anti-squark production and can be comparable to 
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gluino-gluino production depending on the precise squark-gluino mass configuration. 

The outhne of this paper is as follows: In section ^ we review the various tree-level 
contributions to squark-squark production and introduce some notation used throughout 
this paper. Section ^ shows the details to the NLO calculation of 0{a1a) and the strategy 
of the calculation. In Section ^ we list the input parameters used in our numerical analysis 
and show hadronic cross sections and distributions for squark-squark production in proton- 
proton collisions at the LHC. Different SUSY scenarios are considered and a scan over 
squark and gluino masses is performed. Analytic formulas for the tree-level cross section 
and expressions for the cross section in the soft and collinear singular region, as well as the 
Feynman diagrams for the NLO calculation are collected in the Appendix. 



2. Classification of processes and tree-level cross sections 

We consider the pair production of two squarks or two anti-squarks, 

PP^qaq'p, PP^rjp. q,q' = {u,d,c,s}- (2.1) 

where a, /3 = {L, R} label the chirality of the squarks, neglecting left-right mixing. At 
lowest order in QCD there is only one partonic channel for each process, 

qipi)q'{p2) qa{p?,) q'pipi), q{pi) q'{p2) qliPs) q'^(.P4), (2.2) 

where the initial-state quarks and the final-state squarks have to have the same flavor. We 
thus do not consider the production of top (bottom) squarks due to the vanishing (small) 
density of the corresponding quarks inside the proton. Moreover, b-squark production 
has specific different features and will be discussed separately. The unpolarized cross 
sections for squark-squark and anti-squark-anti-squark production are related by charge- 
conjugation. In the following we will refer to squark-squark production only, while the 
charge conjugated processes are properly taken into account in the numerical results. Since 
the electroweak interaction is sensitive to flavor and chirality, one has to treat processes 
with final-state squarks of different chiralities or of different isospin separately, even in the 
limit of degenerate squark masses. CKM mixing effects are neglected in our discussion. 

In total we distinguish 36 processes, resulting from the various combinations of squarks 
of different flavor or chirality in the final state. They can be classified as follows: 

—production of two squarks of the same flavor, 

V , (2.3a) 

PP — )• UaU/3, dadi3, CaCp, SaSjS, {oif3} = {LL, RR, LR}. 

—production of two squarks belonging to the same SU(2) doublet, 

PP — > Uadi3, CaSp, {a/3} = {LL, RR, LR, RL}. 

—production of two squarks in different SU(2) doublets, 

(2.3c) 

PP — )• UaC/3, UaSjs, daCp, daSp, {a/3} = {LL, RR, LR, RL}. 

The corresponding tree-level diagrams of both QCD and EW origin are listed in Figure |l|. 
QCD diagrams are of 0{as), mediated by gluino exchange. EW diagrams are of 0{a) and 
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Figure 1: Parton-level Feynman diagrams for the three classes of squark-squark production at 
tree-level, where a, /? = {L, R}. The first class describes the production of two squarks of the same 
flavor, the second class that of two squarks of the same isospin doublet (but different flavor) and 
the third class refers to the production of two squarks belonging to different isospin doublets. In 
the third class, the subprocess in brackets cannot interfere with other diagrams due to different 
initial state particles. In all three classes, the final-state squarks are of the same generation as the 
initial-state quarks. 



mediated by neutralino or chargino exchange. Quarks and squarks are of the same flavor, 
also in the EW diagrams. The only exception is given by the two pure-EW chargino- 
mediated subprocesses ud — )■ diCL and cd — )• ulsl belonging to the third class, which 
contribute to dici and ulsl final states, respectively. Note that only t- and u-channel 
diagrams are present, but no s-channel diagrams. 

The appearance of both t- and n-channel diagrams for chirality-diagonal QaQa produc- 
tion gives rise to nonzero interferences between QCD and EW diagrams already at tree- 
level.^ The full tree-level contributions to the cross section are thus given by the O(a^) 
Born contribution and the ©(osCe-l-a^) EW contributions. Photon-induced squark-squark 
production is not possible at lowest order from charge and color conservation. 

To keep track of the corresponding order in perturbation theory of the various contri- 
butions, we introduce the notation da"'' ^ [Ai""' in order to refer to the cross section [matrix 
element] at a given order 0{a'^a^) in the strong and electroweak couplings, respectively. 
Results are given in terms of the Mandelstam variables, defined as usual, 

s = {pi + p2)^ i = {pi - psf, u = {pi- Pif. (2.4) 

The differential partonic cross section for a given subprocess qq' — t- qa(i'p at LO can 
thus be written as 



^In the non-diagonal case, qLq'n production, the interference contributions vanish as a consequence of 
the trivial squark mixing matrices in the limit of no L-R mixing, see also the discussion in Appendix 
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in terms of the squared lowest-order matrix element, Ai^'^, averaged (summed) over initial 
(final) state spin and color. Similarly, the pure EW differential cross section of 0{a^) and 
the EW-QCD (D{asa) interference contribution are given by 

da''\s) = Y.\M''^f , (2.6a) 

da''Hs)=T^Re{{M^''yM^'^}^, (2.6b) 

where Ai^'^ denotes the EW tree-level amplitude. Explicit expressions for the squared 
matrix elements are given in [18,24]. For completeness, we include a list of all tree-level 
differential cross sections in Appendix correcting in particular a wrong color factor 
of [24]. 

The hadronic cross sections are obtained from the partonic cross sections by convolu- 
tion with the respective parton luminosity function. At O(a^a^), it is given by 

da'''\S)= [\t ^da'^'^s), (2.7) 
Jto d-T 

with ^ = rTvr^[^^^G'^^) /.'(^'^^) + /.^(^'^^)/.' G'^^ 

Here tq = {mg^ + rriqi^Y / S is the production threshold, determined by the masses of the 
two squarks mq^ and rrifji^. The parton distribution functions (PDFs) fq{x,iJ,p) give the 
probability to find a parton q with momentum fraction x inside hadron A at a factorization 
scale ^F- At the LHC, both hadrons A,B are protons P. S and s = tS are the squared 
center-of-mass (c. m. ) energies of the hadronic and partonic processes, respectively. 

3. Virtual and real corrections of C>(q:^q;) 

At 0{a1a), squark-squark production gets contributions from virtual corrections, real 
photon- and gluon emission, as well as real quark radiation. Ultraviolet (UV) as well as 
infrared (IR) and collinear singularities arise in the one-loop diagrams, see Section The 
IR singularities cancel in sufficiently inclusive observables once virtual and real photon 



and gluon bremsstrahlung corrections are added (see Section 3.2). Remaining collinear 
singularities are universal and can be absorbed by redefining the PDFs, as described in 
Section 

Diagrams and corresponding amplitudes are generated using FeynArts [31,32]. The 
algebraic simplifications and numerical evaluation is done with help of FormCalc and Loop- 
Tools [32,33]. IR and collinear singularities are regularized by means of mass regularization, 
i.e. we introduce a fictitious mass for the photon and the gluon. Quarks are treated as 
massless, except where their masses are needed as regulators. 

3.1 Virtual corrections 

The virtual contributions are given by the interference of tree-level and one-loop diagrams. 
In practice three types of interferences occur at O(a^a), as schematically depicted in 
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EW 




Figure 2: Sample of Fcynman diagrams to illustrate the virtual contributions at 0{a1a). Three 
gauge invariant subsets of interferenees occur at this order. The label of perturbative order is 
attached to each diagram. EW refers to electrovifeakly interacting particles and QCD to strongly 
interacting particles in the loop insertions. The full sets of diagrams are shown in Figs. |l^, 
and ^ 



Figure ^. All three interference terms yield non- vanishing contributions to the cross section. 
For each subprocess, the partonic cross section can be written as 



The first line corresponds to (a) and (6) of Figure y and is given by the interference of 
A^^'*^ with A4^'^. The amplitude A4^'^ is split into two parts, A^J^^^ and -^(qcd)' 
first arising from tree- level QCD diagrams with EW insertions (Figure |2|a, right), and the 
latter from tree-level EW diagrams with QCD insertions (Figure ||b, right). The second 
line in Eq. (3.1), corresponding to Figure |2|c, is given by the interference of Ai'^' ^ with the 



pure-QCD one loop amplitude Ai'^'^. Care has to be taken with diagrams containing a 
four-squark vertex. This vertex includes the electroweak as well as the strong coupling and 
the appropriate part has to be selected in each interference contribution to match the right 
order, as indicated in Figure 

The full set of virtual corrections is UV finite after renormalization of the theory and 
the inclusion of the proper set of one-loop counterterms. The renormalization for squark- 
squark production proceeds in close analogy to that for squark-anti-squark production 
described in [27] and is sketched here only briefiy. Each of the three interference subsets is 
gauge-independent by itself and can be renormalized separately. 

In the first group, shown in Figure §a, UV singularities only arise from gluino-mediated 
amplitudes with weak insertions (A^^^^p. We include the diagrams with counterterms for 



the qgqa vertex, see Figure 13, and evaluate the renormalization constants at 0{a). At this 



order in the perturbative expansion we need to renormalize quark and squark fields, while 
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the renormalization of gluino and strong coupling is not required. The regularization of the 
divergent amphtudes in this sector is done in dimensional reduction, and renormalization 
of quarks and squarks is performed in the on-shell scheme. 

In the second case, Figure ^3, neutralino- or chargino- mediated amplitudes with strong 
insertions (-^(qcd)) are considered. To obtain a UV-finite result, one needs to include di- 
agrams containing counterterms for the qqx'^ vertex and, if arising, for the qq'x^ vertex. 



see Figure 14. The renormalization constants have to be evaluated at 0{as) and no renor- 
malization of the neutralino or chargino is required. Since the gluino does not enter this 
subset of one-loop amplitudes, it is thus sufficient to renormalize the quark and squark 
sector. As before, the divergent amplitudes are regularized in dimensional reduction and 
on-shell conditions are imposed to fix the (s)quark renormalization constants. 

The third subset, Figure |^, refers to pure-QCD one-loop amplitudes, i. e. gluino- 
mediated diagrams with strong insertions (A^^'"). In this case one has to renormalize the 
quark and squark sector as well as the gluino and the strong Yukawa coupling gs, which 
appears in the qqg vertex. The renormalization constants in the corresponding amplitudes, 



see Figure 15, have to be evaluated at 0{as). The strong scalar coupling gs is related to 
the strong coupling g^ via super symmetry. To match the definition of the strong coupling 
constant used in the extraction of the PDFs, gs has to be given in the MS scheme with the 
contributions from heavy particles subtracted in the running of Og ■ We thus regularize this 
part of the virtual corrections using dimensional regularization. Quarks and squarks are 
renormalized on-shell again, in the strong sector the MS scheme is applied. Dimensional 
regularization however induces a finite difference between gs and gs at the one-loop level 
and violates the super symmetric relation between the two couplings [34]. We add the well- 
known finite shift in the definition of the renormalization constant for gs in order to restore 
SUSY in the physical amplitudes [18]. 

All counterterms and renormalization constants are explicitly given in [27], Appendix 
B, and need not to be repeated here. 

The virtual corrections also exhibit photonic and gluonic mass singularities of infrared 
(IR) and collinear origin. In A^^^^^ mass singularities arise if two external particles 
exchange a low-energetic massless photon while collinear singularities appear if one of the 
massless initial-state quarks splits collinearly into a quark and a photon. In order to 
obtain an IR finite result, real photon radiation at ©(a^a) has to be added. In contrast 
in A1(QCJ3)) massless gluons running in the loops give rise to mass singularities in the soft 
and collinear limit. Similarly, the diagrams contributing to Ai'^'^ suffer from gluonic IR 
and collinear singularities. Hence we have to include real gluon bremsstrahlung at 0{ala) 
in order to cancel the IR singularities. We regularize the photonic singularities by means 
of mass regularization. Owing to the photon-like appearance of the gluon in the respective 
diagrams, it is also possible to regularize these IR singularities by a fictitious gluon mass. 

3.2 Real corrections 

Three independent bremsstrahlung processes contribute at 0{ala), as depicted in Figure |3|. 
Real photon and real gluon radiation processes have to be combined with the corresponding 
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photon bremsstrahlung 



gluon bremsstrahlung 



real quark radiation 




Figure 3: Sample of Feynman diagrams for the three subsets of real emission contributions at 
0{ala). The order in the perturbative expansion is specified for each diagram. The full sets of 
diagrams are given in the Appendix, cf. Figures |l6|, |l^, and |l^. 



subset of virtual corrections to obtain an IR finite result. Also real quark radiation gives 
nonzero contributions from the interference of QCD and EW mediated diagrams and has 
to be included in the cross section at 0{a1a). 

3.2.1 Real photon emission 

The real photon emission at 0{a1a), 

q{pi)q{p2) qa{p'i)qfj{PA)i{k), (3.2) 

is given by the squared matrix element of a QCD tree-level diagram with an external photon 
attached (see Figure ^, top and Figure ^ for the full set of diagrams) . The integration over 
the photon phase space is IR divergent in the soft-photon region, i. e. for A;'' — )• 0. Further 
singularities arise in the collinear region if • — )• for i = {1,2}. We use phase space 
slicing and apply a cut on the photon energy, > Ss\^/2, and on the angle 6 between 
the photon and incoming partons, | cos{9)\ < 1 — 6$, to split off the singular regions. In 
the hard, non-collinear region the integration is convergent and is performed numerically. 
The cross sections in the soft- and collinear region can be approximated analytically and 



are given in Appendix B.l 



By combining the real photon emission with the virtual EW-type corrections (Fig- 
ure ^), the soft singularities cancel. Remaining initial-state collinear singularities are 
universal and have to be absorbed via factorization in the PDFs, see section |3.3| . In the 
following, we will refer to this UV-, IR- and collinear-finite combination as the EW-type 
corrections. 

3.2.2 Real gluon emission 

Real gluon bremsstrahlung at O(a^a) proceeds via the partonic process 

q{pi)Q{p2) qa{p'i)qp{p?,)g{k). (3.3) 
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It is given by the interference term of a QCD- and an EW tree-level diagram, both with 
an external gluon attached on (Figure ^, center and Figure 17). In the considered pro- 
cesses the gluon is Abelian like and we can treat soft and collinear singularities by mass 
regularization in close analogy to the photonic case. However the eikonal current has to be 
modified in order to take color correlations into account. Different to real photon emission, 
collinear singularities only arise for diagonal QaQa, ulcIl, and clsl production. This can 
be seen by noticing that in the collinear cone the cross section becomes proportional to 
the corresponding 2 — )■ 2 process, which in this case would be squark-squark production at 
0{asa); i.e. the interference of tree-level QCD and EW diagrams must be non-vanishing. 
Explicit expressions for the cross sections in the soft- and collinear regions are given in 
Appendix B.2 . 

By combining real gluon emission and the two virtual QCD-type corrections (Fig- 
ure P3,c), the IR singularities cancel. Remaining collinear singularities are again absorbed 
into the PDFs, as described in Section |3.3| . In the following we will refer to this UV-, IR- 
and collinear-finite combination as the QCD-type corrections. 



3.2.3 Real quark emission 

Finally, also real quark radiation contributes at ©(a^a), 

9{pi)qiP2) qa{p^)q'f}{pA)q{k), ^^^^ 

and iiq^q g{pi) q{p2) ^ qaiPs) q'piPi) q{k), 
via the interference of a QCD-type diagram with an EW-type diagram, as shown in Fig- 



ure y, bottom (see Figure 18 for the complete listing of diagrams). This process can 
be regarded as completely independent to the virtual corrections, since it is IR finite by 
itself. However it has to be taken into account in a consistent analysis of electroweak cor- 
rections up to 0(0^0). Initial-state collinear singularities are present for processes with 
non-vanishing tree-level interferences. The cross section in the collinear region is given in 



Appendix B.3. As before, these singularities are absorbed via factorization into the PDFs, 
see Section |3.3| . 

Different to photon and gluon bremsstrahlung, the internal gluino, neutralino or charg- 
ino can go on-shell in specific SUSY scenarios, if heavier than one of the external squarks. In 
these cases, we include a Breit-Wigner width for the resonant particle in the corresponding 
propagators to regularize the poles. Note that physical resonances do not occur. This is 
different to the case of real quark radiation in e. g. gluino-squark production processes [28], 
where internal squarks can go on-shell in both the EW- and the QCD-mediated diagrams. 



3.3 Factorization of initial-state collinear singularities 

The remaining collinear singularities have to be absorbed by redefining the PDFs. At 
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0{a1a) this can be achieved by the replacement [35,36] 

/ aei + agCp 1 ae^ 

fq{x,^F) fq{x,flF) 1 Ky+s " T fvA 

\ TT 4 vr 



z -^^Vz'""; \ 27r ■ 27r 

where denotes the electric charge of quark q, Cp = 4/3, and 
Kv+s = 1 - In (5s - In^ 6s + {^n 5., + - ] In 

The factorization-scheme dependent functions are 

27r2 n 
= 9 + — + 31n5s-21n2 5„ 



4 y V 



Uz) = Pqq{z) In (1-^) - + 2Z + 3, 



(3.6) 



(3.7) 



with the splitting functions 

P,,i^) = ^-^^ P,,{z)=z' + il-z'). (3.8) 

The factorization is done in the MS scheme at NLO QCD and in the physical DIS scheme at 
NLO EW. The replacement of the PDFs in Eq. (2/7) gives further contributions of ©(a^a) 
to the total cross section. The first and second line in Eq. (3^) cancel the remaining 
singularities in the EW-type and QCD-type corrections. The third line in Eq. (|3.5D cancels 
the collinear singularities in the real quark radiation. 



4. Numerical results 



In the following we illustrate the impact of the EW contributions on the production cross 
section. Since we have 36 processes contributing to squark-squark production and the 
same amount for anti-squark-anti-squark production, we present (at least partly) inclusive 
results. We refer to four different combinations of (anti-)squarks in the final state, which 
differ with respect to the chirality of the produced particles: 

• qiQ^ or "LL" refers to the inclusive production of two left-handed squarks and two 
left-handed anti-squarks. 

• qLQR or "LR" refers to the inclusive production of one left-handed and one right- 
handed squark and the charge conjugated process. 
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• qnq'fi or "RR" refers to the inclusive production of two right-handed squarks and two 
right-handed anti-squarks. 

• qq' or "inch" refers to the inclusive production of all (anti-)squarks. It is given by 
the sum of the three cases above, taking all 72 subprocesses of squark-squark and 
anti-squark-anti-squark final states into account. 

We focus here on these chirality-based classes since squarks of different chiralities are, in 
principle, experimentally distinguishable by their decay chains, see e. g. Section 5.1.2 of [37]. 

In the discussion we refer to the following quantities, based on the cross section def- 
initions in Section |2[ The leading order cross section is denoted by a^°'^'^ = a"^'^. The 
tree-level EW and the NLO EW contributions to the cross section are labeled by 

respectively, and Act^^ = Act*''^*' + Aa^^'^ wih be referred to as the EW contri- 
bution. The total sum of the LO cross section with the EW contributions is denoted by 
o-NLO ^ ^Born j_ ^^EW^ Relative EW contributions are defined by 

jtree EW _ ^^trec EW ^^Born ^NLO EW _ ^^NLO EW ^^Born ^EW _ ^^EW ^^Born 

(4.2) 

In distributions 6 denotes the relative EW contribution defined as 5 = (Onlq— C'Born)/C'Born, 
where O is a generic observable and Onlo is the sum of the Born and the EW contribution. 

4.1 Input Parameters 

The Standard Model input parameters are chosen in correspondence with [38,39], 

Mz = 91.1876 GeV, Mw = 80.4247685 GeV, 

a"^ = 137.036, a^(Mz) = 0.119, (4.3) 

mt = 170.9 GeV, m^^ = 4.7 GeV. 

The strong coupling constant has been defined in the MS scheme using the two-loop 
renormalization group equation with five light fiavors. 

For the SUSY parameters, we refer to three benchmark mSUGRA scenarios, the SPSla' 
scenario, the SPS2, and SPS5 scenario [38,40]. The SPSla' scenario can be considered as a 
"typical" mSUGRA scenario. It has been proposed by the SPA convention and should be 
used for comparisons with other calculations. The SPS2 scenario features relatively heavy 
squarks with light charginos and neutralinos and a gluino lighter than the squarks. The 
SPS5 scenario leads to a very light ti with moderate light-flavor squark masses. In each sce- 
nario, the particle spectrum is determined by universal GUT scale parameters, cf. Table |l|, 
which act as boundary conditions for the renormalization group running of the soft-breaking 
parameters. We use the program Softsusy [41] to evolve the soft-breaking parameters down 
to the SUSY scale MsusY- We choose a common SUSY scale MsusY = 1 TeV for all sce- 
narios, in reference to the SPA convention. At MgusY a consistent translation of the squark 
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mo 


mi/2 


Ao 


tan P 


sign(/i) 


SPSla' 


70 GeV 


250 GeV 


-300 GeV 


10 


+ 


SPS2 


1450 GeV 


300 GeV 





9.66 


+ 


SPS5 


150 GeV 


300 GeV 


-1000 GeV 


4.82 


+ 



Table 1: High energy input parameters for the different SUSY scenarios considered. The mass 
parameters toq, mi/2 and Aq arc given at the GUT scale, tan/3 is evaluated at Msusy = ITeV. 





ul 


ur 


dL 


dR 


a 




xt 


SPSla' 


561 


543 


566 


539 


609 


101 


180 


SPS2 


1559 


1554 


1561 


1555 


785 


120 


199 


SPS5 


677 


655 


681 


654 


724 


123 


225 



Table 2: On-shell masses of the squarks, the gluino, and the lightest neutraliuo and chargino 
within the different SUSY scenarios considered. All masses are given in GeV. 

PP -> {qLqL,qRqR,ULdL,CLSL} 

420 

400 

380 

1^360 
b 

340 
320 
300 

0.5 0.6 0.7 0.8 0.9 1 2 

/i/(560GeV) 

Figure 4: Dependence of the hadronic cross section a on the rcnormalization and factorization 
scale /i/i and hf for the SPSla' scenario. Both scales are set to a common value = /i^ = A* 
with /z varied by a factor two around = 560 GeV. Only processes that contribute at 0{asa) 
are considered. The red line shows the EW tree-level contribution to the cross section, while the 
green curve shows the EW contribution up to NLO EW. For the latter case the scale dependence 
is reduced. 

masses into the on-shell scheme is performed. The left-handed down-type squark is treated 
as a dependent mass parameter, fixed by SU(2) invariance. It is set to its corresponding 
on-shell value obtained at one-loop accuracy according to [42]. The on-shell mass param- 
eters for the light flavor squarks together with the masses of the gluino and the lightest 
neutralino/chargino are summarized in Table ^. 

The technical cuts needed for the regularization of soft and collinear singularities are 
set to 6s = 10~'^\/I and 60 = 10~^. We checked numerically that these values are sufficiently 
small to justify the eikonal approximation. For the calculation of hadronic observables we 
use the MRST2004QED parton distribution functions [43]. The factorization and renor- 
malization scales are set to a common value, = fin = fip = rriq with rriq being the average 
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SPSla' 


^Born 


^^tree EW 

0{asa + a^) 


^^NLO EW 


^tree EW 


^NLO EW 


^EW 




1717.6(8) 


378.9(1) 


-74.8(6) 


22.1 % 


-4.4 % 


17.7 % 




1981.9(7) 


31.81(2) 


-1.60(9) 


1.6 % 


-0.1 % 


1.5 % 




1743.8(4) 


2.538(1) 


-70.71(4) 


0.1 % 


-4.1 % 


-3.9% 


qq' 


5443(1) 


413.3(1) 


-147.1(6) 


7.6 % 


-2.7 % 


4.9 % 



Table 3: Hadronic cross sections for squark-squark production at the LHC within the SPSla' 
scenario. Shown are the LO cross section, the tree-level EW as well as NLO EW contributions and 
the relative corrections as defined in the text. Anti-particles are included. The numbers in brackets 
refer to the integration uncertainty in the last digit. All cross sections are given in fcmtobarn (fb). 



SPS2 


^Born 


^^tree EW 

0{asa + a^) 


^^NLO EW 


^tree EW 


^NLO EW 


^EW 




7.359(1) 


1.0326(2) 


-0.5776(7) 


14.0 % 


-7.8 % 


6.2 % 




7.529(1) 


0.1005(1) 


-0.0052(1) 


1.3 % 


-0.1 % 


1.3 % 


QLq'u 


14.651(1) 


0.0136(1) 


-0.8676(2) 


0.1 % 


-5.9 % 


-5.8 % 


qq' 


29.539(2) 


1.1468(2) 


-1.4506(7) 


3.9 % 


-4.9 % 


-1.0 % 


Table 4: Same as Table ^ but for the SPS2 scenario. 


SPSS 


^Born 

0(«^) 


^^tree EW 


^^NLO EW 


^tree EW 


^NLO EW 


^EW 




774.7(1) 


185.71(4) 


-35.9(1) 


24.0 % 


-4.6 % 


19.3 % 




888.0(1) 


16.332(5) 


-0.69(2) 


1.8 % 


-0.1 % 


1.8 % 


QLq'u 


758.00(9) 


1.1559(3) 


-33.68(1) 


0.2 % 


-4.4 % 


-4.3 % 


qq' 


2420.7(3) 


203.20(4) 


-70.3(1) 


8.4 % 


-2.9 % 


5.5 % 



Table 5: Same as Table | but for the SPS5 scenario. 



mass of the light-flavor squarks. Figure |^ shows the scale dependence for the SPSla' sce- 
nario. Here we are inclusive with respect to processes that contribute at 0{asa). We 
compare the tree-level EW contribution with the cross section including the EW NLO 
contributions at O(a^a), i.e. one order higher in the strong coupling. One finds that the 
scale dependence is considerably reduced. 

4.2 Total hadronic cross sections 

Tables |3|-|5| give the results for the hadronic cross sections for squark-squark production 
at the LHC within the SPSla', SPS2, and SPS5 scenario, respectively. We refer to the 
production of squarks of different chiralities separately. Renormalization and factorization 
scales are set to /i = 560 GeV (SPSla'), // = 1560 GeV (SPS2) and /i = 666 GeV (SPS5). 
The Born cross section is QCD mediated and does not depend on the chirality of the 
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produced squarks. Indeed the cross sections for the diagonal production of two squarks, LL 
or RR, become equal for degenerate masses. The Born cross section for non-diagonal LR 
production, however, is different in general since the final-state particles are distinguishable. 
In the SPSla' and SPS5 scenario it happens to be of similar size as the LL and RR 
production cross sections. In the SPS2 scenario, however, the LR cross section is enhanced 
and accounts for 50% of the total cross section for inclusive squark-squark production. As 
we will see below, cf. Figure |^(b), the relative yield of the LR production cross section 
is determined by the ratio of squark and gluino masses, becoming more important if the 
exchanged gluino is lighter than the final-state squarks. 

Due to the nature of the electroweak interaction, the EW cross section contributions 
depend strongly on the chirality of the produced squarks. The impact of tree-level EW 
contributions is largest in case of LL production (15 — 25% in the considered scenarios) and 
around one order of magnitude smaller for RR production. It is even further suppressed 
for LR production, where the 0{asa) tree- level interference contributions are completely 
absent. The situation is different for the NLO EW contributions. These are equally 
important in the case of LL and LR production (reducing the LO prediction by about 4—8% 
in the considered scenarios), but negligible for RR production. In all three cases the NLO 
EW corrections are negative and partially compensate the EW tree- level contributions. 
Summing over all processes, the EW contributions to the total cross section for inclusive 
squark-squark production decrease from about 8% to about 5% in SPSla' and SPSS after 
the inclusion of NLO EW corrections. In the SPS2 scenario, where LR production is the 
dominant production mechanism, the NLO EW corrections even over compensate the EW 
tree-level contributions and the result turns negative. 

In order to further investigate the dependence of the EW contributions on squark and 
gluino masses, we perform a parameter scan on those quantities. The independent squark 
masses are chosen equal to a common value rriq, with the dependent fourth squark mass set 
to its corresponding on-shell value. All other parameters are fixed to their SPSla' values. 
The renormalization and factorization scale are identified with the common squark mass 
AiiJ = AiF = ruq. 

To start with, we show in Figure ^ the Born cross section (left panel) and the EW cross 
section contributions (right panel) for inclusive squark-squark production as a function of 
the common squark mass and for different values of the gluino mass. Both the Born and 
the EW contribution strongly decrease for growing squark masses. While the Born cross 
section is quite sensitive to the gluino mass for low squark masses, the EW contribution is 
almost independent in this regime. For high squark masses the behavior is vice versa. 

Figure I shows the K factor, defined as K = a^^'^/a^°"^, for the same parameter 
range as considered in Figure |5|. The three different combinations of chiralities of the 
final-state squarks, as well as the inclusive case are considered separately. The K factor 
is largest for two left-handed squarks in the final state. Here, the EW contributions alter 
the LO cross section prediction between 10 — 50%, being most important in case of light 
squarks and a heavy gluino. The EW contribution is enhanced by the large tree-level EW 
contribution. In the case of RR production, the EW contributions are below 3% in most 
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Figure 5: Hadronic Born cross section (left) and EW contributions (riglit) for inclusive squark- 
squark production as a function of a common squark mass mq. Different gluino masses rrig are 
considered, all other parameters are set to their SPSla' values. 

parts of the parameter space. For LR production the EW contributions are mainly given by 
the NLO EW corrections, leading to a K factor smaller than unity. The LO cross section 
is reduced by —3% to —5%, most strongly in a scenario with heavy squarks and a heavy 
gluino. Altogether one finds for inclusive squark-squark production EW contributions that 
range from 9% for rrig = 400 GeV up to 22% for rrig = 1000 GeV for light squarks. For 
heavy squarks, the EW contributions are only at the percent level due to the interplay of 
positive EW corrections in the LL and RR case and negative EW corrections in the LR 
case, suppressing the EW contributions by one order in magnitude. 

We can understand the smallness of the EW contribution for high squark and low 
gluino masses by having a closer look at the interplay of the tree-level EW and NLO EW 
contributions. The ratio Ao"^^*-* EWy^^tree EW £qj, inclusive squark-squark production is 
shown in Figure 0(a). As one can see, the NLO EW corrections become more important 
for larger ratios raq/rrig and reach the same size as the tree- level EW contributions for 
about mq/rrig > 1.5, depending on the precise value of the gluino mass. This is due 
to the fact that the LR contribution becomes more relevant for increasing rrig/mg, see 
Figure |^(b). Owing to the suppressed tree-level contributions, the EW contributions to 
LR production are negative and partially compensate the positive yield from LL and RR 
production. Figure |^ (a) also confirms our observation from the SPS2 scenario that the 
NLO EW corrections compensate the tree-level EW contributions in the inclusive cross 
section, cf. Table which seems to be a generic feature in scenarios with the squark 
heavier than the gluino. 

4.3 Differential distributions 

Here, we illustrate the results for the SPSla' scenario. In Figures I, I, and ^ we consider 
the differential distributions of the EW contributions with respect to various kinematical 
variables. In the left panels, the tree-level EW contributions and the three gauge-invariant 
subsets of NLO EW contributions (EW-type corrections, QCD-type corrections, real quark 
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Figure 6: K factor defined a.s K ^ a^^'^ / <t^°"\ as a function of a common squark; mass. Different 
gluino masses rrig are considered, all other parameters are set to their SPSla' values. The labels 
"incl." and "LL" , "RR" , "LR" refer to inclusive squark-squark production and chirality-grouped 
subprocesses as explained in the text. 



radiation), as well as the summed EW contributions are shown. In the right panels, the 
impact of the EW contributions relative to the Born cross section, 5, is given. 

Figure |8| refers to the distribution with respect to the transverse momentum pT of 
the squark with highest pT- The tree- level EW contributions are always positive with a 
maximum at about 250 GeV and dominate the sum over a wide range of the phase-space 
for LL, RR, and inclusive squark-squark production. Again, they are suppressed for LR 
production. The interplay of the NLO EW contributions is more complicated. For all 
processes, the real quark radiation is small and mostly negative. For LL production, large 
cancellations among the EW- and QCD-type corrections occur. As a result, the relative 
yield is dominated from the tree-level contributions in the small-pT region where it is large 
and positive (up to 25%). For higher values of pT, the relative corrections turn negative and 
grow up to —10%. In case of RR production the EW-type corrections are suppressed from 
the chirality and the QCD-type corrections are more important. However the relative EW 
contributions in total do not exceed a few percent. Finally in the LR case, the QCD-type 
corrections are negligible since they are related to QCD-EW interferences. The dominant 
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contribution arises here from the EW-type corrections. The relative contributions are 
always negative, between —2% for small values of px and up to —10% in the high-px 
region. It is important to note that even though the relative NLO EW contributions to 
the integrated cross section are comparable for LL and LR production, cf. Table ^ they 
originate from distinct sources and the differential distributions differ strongly. 

In Figure ^ the distributions with respect to the invariant mass of the squark pair, 
= ips + Pi)'^, are shown. The interplay of the various subsets of EW contributions is 
similar as for the px distributions. For LL and RR production, the EW tree-level contri- 
butions are dominant and peak at around M[^^ ~ 1200 GeV. They shift the relative EW 
corrections to positive values, up to 30% in the low-Mjnv region for LL production. In case 
of non-diagonal LR production, where the EW tree-level contributions are suppressed, the 
relative corrections are negative and grow up to —5% for the intermediate and high-energy 
region. Finally, in the inclusive case, we find a strong energy dependence of the relative 
EW corrections, ranging from +10% for Mjnv ^ 1200 GeV to -5% for Mi^v > 3500 GeV. 

In Figure ^ we present the pseudo-rapidity distributions, where always the squark with 
the higher absolute value of the pseudo-rapidity r] (in the laboratory frame) is considered. 
All EW contributions are vanishing in the central 77 ~ region. The characteristics of 
the rapidity gap in the distributions depend on the precise quantity considered and is 
enhanced by our choice of referring to the larger 77. The tree-level contributions peak at 
around |7/| = 2 and dominate the total result, if present. The distributions for EW-type 
and QCD-type NLO corrections differ in sign and shape from each other, leading to large 
cancellations over a wide phase-space range. In total, the EW contributions alter the LO 
distributions by up to 20 — 40% for LL and up to 10 — 20% for inclusive squark-squark 
production in the strong forward region for \rj\ > 2. 
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Figure 8: Transverse momentum distributions for squark-squark production at tlie LHC within 
the SPSla' scenario. Shown are the tree-level and NLO EW cross section contributions (left) and 
the impact of EW contributions relative to the QCD Born cross section (right) for inclusive qq' 
production (top), production of two left-handed squarks Qlq'l (second), production of two right- 
handed squarks qnq'fj (third), and non-diagonal qLQR production (bottom). Charge conjugated 
processes are included. 
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Figure 9: Invariant mass distributions for squark-squark production at the LHC within the SPSla' 
scenario. Shown are the tree-level and NLO EW cross section contributions (left) and the impact of 
EW contributions relative to the QCD Born cross section (right) for inclusive qq' production (top), 
production of two left-handed squarks qLtj'^ (second), production of two right-handed squarks qRq'fj 
(third), and non-diagonal qLq^ production (bottom). Charge conjugated processes are included. 
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Figure 10: Pseudo-rapidity distributions for squark-squark production at the LHC within the 
SPSla' scenario. Shown are the tree-level and NLO EW cross section contributions (left) and 
the impact of EW contributions relative to the QCD Born cross section (right) for inclusive qq' 
production (top), production of two left-handed squarks qh^L (second), production of two right- 
handed squarks qnq'j^ (third), and non-diagonal qL<i'ft production (bottom). Charge conjugated 
processes are included. 
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Figure 11: (a) Transverse momentum distribution of the hardest squark for ULd-L production 
within the SPSla' scenario. Strong cancellations occur between the different contributions at NLO 
EW. (b) Relative NLO EWcontribution, defined as the ratio of dcr^LO ew/^^^ .^^^ d(T^°"'/dpT. 



Up to now, our discussion has only been for inclusive combinations of final-state squarks 
for given chiralities. To get further insight on the cross section, we show in Table ^ the 
cross section divided up into the various subprocesses for squark-squark production within 
the SPSla' scenario. Again, anti-particles are included. Owing to the degenerate masses 
of first- and second-generation squarks, we do not distinguish between final states that 
result from exchanging both squarks with their first or second generation counterpart, 
i.e. ulul production also includes clcl production, etc.. This reduces the number of 
distinct subprocesses from 36 down to 22. The contributions to Ao"'''''"' are always 
positive and are largest for uidi production due to the interference of gluino and chargino 
exchange diagrams and constitute 57% of the inclusive tree-level EW contribution, see also 
Table ^. One even finds that the inclusive tree-level EW contribution is given to 98% 
by only five processes, namely uiui, urur, didL, uidi and ulsl- The contributions to 
^^NLO EW j-Qostly negative, reducing the importance of EW contributions. In contrast 
to the tree-level EW case, many processes contribute with a significant amount to the 
inclusive NLO EW contribution of the cross section. Especially for processes with squarks 
of different generations, Aa^^ is mostly dominated by NLO EW contributions. The size 
of the NLO EW contributions is often reduced due to the interplay of QCD-type and EW- 



type corrections as shown in Figure 11 in the case of ULdi production. The different types 
of NLO EW corrections partially cancel. Furthermore, the sum contains corrections of 
positive and negative sign, leading to an integrated result Ao"'^'"*-* that is considerably 
smaller than the corrections affecting the LO result in various phase-space regions. 

5. Conclusions 

We have studied the hadronic production of two squarks or two anti-squarks within the 
MSSM, including tree-level EW and NLO EW cross section contributions of 0{as a + a^) 
and ©(a^ct), respectively. In contrast to pure QCD computations, one has to treat pro- 
cesses with squarks of different fiavor or chirality separately. 
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Table 6: Hadronic cross section for squark-squark productfon at the LHC within the SPSla' 
scenario. Charge conjugated processes are included, uu final states include cc, dd include ss, and 
us includes cd. All cross sections are given in femtobarn (fb). 



At 0{a1a), numerous QCD-EW interference terms occur. Virtual corrections arise 
from the interference contributions of tree-level QCD amplitudes and mixed EW-QCD 
one-loop diagrams, as well as from the interference of tree-level EW and pure-QCD one- 
loop amplitudes. Bremsstrahlung corrections comprise real photon, real gluon, and real 
quark radiation processes. 

We have performed a detailed numerical analysis for two left-handed squarks (gLQi), 
two right-handed squarks {qRqR), and one left-handed and one right-handed squark (qLQR), 
as well as for inclusive squark-squark {qq) production at the LHC. The tree- level EW con- 
tributions are largest for qLqL production, where they are enhanced by ©(a^a) interference 
terms and can easily reach the 20% level. The interference contributions are suppressed for 
qRqR production from the chiral couplings and vanish for qiqR. production in the limit of 
no L-R mixing. At NLO, partial cancellations occur among the various EW contributions. 
As a result, the integrated cross section is reduced by the NLO EW contributions by a few 
percent for most subprocesses. The interplay of tree-level and NLO EW contributions is 
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not universal and depends sensitively on the ratio of squark and gluino masses as well as 
on their absolute values. The full EW contributions affect the integrated cross section for 
inclusive squark-squark production at the percent level (about 5% in SPSla' and SPSS, 
— 1% in SPS2). In the distributions, however, the EW contributions range from —10% to 
20% and even larger values for exclusive (JlQl production. 
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A. Tree-level cross sections 

In this Appendix, we give explicit expressions for the tree-level differential cross sections 



Eq. ( ^.SD -Eq. (|2.6D for all squark-squark production subprocesses. We closely follow [24] 
and express the color- and spin-averaged squared t- and n-channel matrix elements and 
their interference in terms of the following functions. 
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Here, ^1^2 £ {g, X^, X^} denote the exchanged particle in the t- or u-channel, respectively, 
and C|$ 0^^} are color matrices which are summarized in Table |^. Note that the color 
factors Cii-(x,x) differ from [24]. 
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Table 7: Color factors in Eqs. |A.1|-|A.3|, with x denoting any neutralino Xk chargino Xfc • 
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a',b' 

^-channel: 2 



a, b 




Figure 12: Notations for the couplings in the tree-level formulas in Appendix ^ Labels a,c refer 
to the couplings to left-handed squarks, b and d to that of right-handed squarks. Couplings at the 
upper and lower vertex, respectively, are denoted by distinct labels for convenience. 



^15 Qa 


a, a', c, c' 


b, b',d,d' 














Xki< 







X^da 


-^U*,,6La 







-V^igs 5 La 




9, da 


-V2igs 5 La 


^igs SRa 



Table 8: Coupling constants a^, ~ q t ■ ■ exchange particle and produced light- flavor 
squark g^, following the conventions of [33]. L-R mixing of the squark mass eigenstates is neglected. 
N, U, V are the unitary matrices diagonalizing the neutralino and chargino mass matrix, respectively 
and s,,, = sin and c,„ = cos 6,,, . 



The coupling constants are collected in the abbreviations A, B, C, D, F, as given by 



Si.Iq i2,qc ii,q'i3 i2,q'fi Ci,-?^ ^2,9^ 51:9^ 52,ga' 



D{ii,i2,qa,q[s) 
F{ii,i2,qa,qB) 



O-e n 0*? - <^'? -/ a?* -/ + b( p; - b'- 6'r* , 

a.'jQ 6,<?a 6. 9^ SI'"?" 6:9a 6:9^ Ca-Qa 

Cl:9;3 6,9^ 6:90 6:9a 6:9a 6:9a 6.93 6:9^ 

-/ Cr -/ Cr ~ Cz - ~\~ ~i dz -., dz - dz ~ . 

6.9,3 6.9J3 6.9a 6.9a 6.9,3 6:9^ 6:9a 6.9a 

a; - c| d- c'-* - -I- 6; - (i| 6'- _ , 

6.9a ^2:9^ 6:9^ 6:9a 6.9a ^2:9^ 6.93 6.9a' 



(A.4a) 



where the notation refers to the labels as listed in Figure 12. Finally the explicit coupling 
constants a? - , 6; = , . . . are given in Table 

si. 9a si. 9a L-l 

For the differential cross sections, we refer to the three classes of subprocesses intro- 
duced in Eq. ( |2.3| ). a, /3 = {L, R} label the chirality of the squarks, /c, I label the four (two) 
mass eigenstates of neutralinos (charginos). 

• PP — )■ (two squarks of the same flavor) 

The partonic process for this class of processes is qq — ?• qaq^^ i- e. all quarks and 
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squarks are of the same flavor. The differential cross sections at ©(a^), 0(0^), ©(a^a) 
read, according to the notation introduced in Section |2|, 



dt 



k.l=l 



(A.5) 

dt 



+ 2Re{*(x2,X?,9a,g/3)} 

4 ^£ 

dai'i = 5]2Re{M'(5,X?,g~a,9»} 

fc=i 

As can be seen from the couplings in Table ^ the interference terms ^ and thus in 
particular the interference contribution da^' ^ are only present for diagonal squark- 
squark production (i. e. a = (3). This is a result from the absence of L-R mixing for 
the light-flavor squarks. 

PP — )• QciQja (two squarks of different flavor in the same SU(2) doublet) 

The only contributing partonic process is qq' — t- qaq'js, with being the SU(2) partner 
of q. The tree-level contributions to the cross section read: 

da'^'^ =<^>{g,g,qa,q'R 



da0'2=| ^ Hxlx'^,qa,q'p)+Y. ^ixlxf ,qa,q'p) 



k,l=l kd=l 



k=l 1=1 

2 



IGvrs^ ' 



da^^^ = Y2Re{^ig,xt,qo^,q'^)} 



fc=i 

In this case, the interference terms are related to chargino-mediated diagrams and 
thus the interference contribution da^' ^ is only non-zero for the production of two 
left-handed squarks (i. e. a = /3 = L). 

PP QoL^'i^ (two squarks in different doublets) 

This class describes the production of two squarks of different flavor and of different 
generation, arising from the partonic process qq' — >• qaq'p, with q / q' . The tree-level 
cross sections read as follows, 

da''^=^{g,~g,qc..q'p)-^, (A.7a) 



- 25 - 



da 



kA=l J 



k,l=l 

di (A.7b) 



da^'^ =0. 



Here, two additional chargino-mediated partonic processes {us — )• d^CL and cd — )• 
ulsl) can give an O(a^) contribution. The 0{asa) interference contribution vanishes 
for this class of processes. 



B. Bremsstrahlung integrals 

Here we give the soft and collinear singular integrals appearing in the phase space integra- 
tion of the real radiation processes, with either a photon, gluon or (anti-)quark radiated. 

B.l Soft and collinear photon bremsstrahlung 

The cross section for real photon emission factorizes in the soft limit from the original cross 
section without photon emission, 



da 



2, 1 

12-^347 



soft 



a 
2^ 



{ Yl eiCjaiajIij^ 



da 



2,0 

12-5-34' 



(B.l) 



i,j=l;i<j 



where now the considered process is given as a subscript for clarification and particles are 
labeled by i = 1 . . .4 according to the definition of momenta pi in Eqs. ( p. 21) and (p.2|). 
Ci is the charge of the ith particle and fij = ±1 depending on whether the particle is 
incoming or outgoing, respectively. The phase-space integrals lij = Iji are universal and 
well-known [44]. In the limit of vanishing initial-state masses (and two massive final-state 
particles), they are given by 



In 



(M^Ui,/'"- 



x..^i..(l(M)!,, 



Sl2 
1 . /I 



In , 

A VI + A 



for i = {1,2}, 
for i = {3,4}, 



12 



i=l,2 



In 



S12V f4{AE) 



In 



S12 
2 



; inM - 



m 



"3 



2n 



J. 



■34 



T- ■ 



1 

V34 



In 



In 



i=3,4 

mfm?- 



In 



1 + A 

1-A 

/4(AE)2 



In^ 



V A2 



1 



S12 



- InM ^ - - In 



2Li2 1 



(1 + /3, 



2 

2Li2 I 1 



rri: 



1-A 



vr 
"3 



1-A 
2 



(B.2) 



4^(1-/3,) 



for 



i = {l,2}, 
j = {3,4}. 
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with Sij = 2pi ■ pj, l3i = Ipil/p'i, Vij = y 1 — 4mimj/sfj, AE = a/Ii25s/2. A is the fictitious 
photon mass. 

The partonic cross section in the cohinear region can be written in terms of a convo- 
lution integral 



= ^- ^ / dz Kqq{z,s) dal'^-^^^izs) , 

coll ivr J^^ 



(B.3) 



where zq and k„„ are given by 



zq = (m? + ml)/s, 



, .^ 1 + , ( sbe 

Kqq{z,s) = In 



1 — z V 2mq J 1 — z 



2z 



(B.4) 



The upper integration bound is lowered by 6s to avoid double counting of the soft regime. 
B.2 Soft and collinear gluon bremsstrahlung 

The bremsstrahlung integrals for soft gluon emission are the same as those for the photonic 
case, but in addition one has to take the color correlation of the amplitudes into account. 
We use the following notation to keep track of the color factors. 

Let I ci, . . . , Cm) denote a complete color basis. The colored matrix element with m external 
particles Pi carrying momentum pi and color a is then given by 



. /JJ Ci...Cm 

-'^'12-H>3...m 



Cl , . . . , Cji 



(B.5) 



For the color structure in the case of real gluon emission it is convenient to associate a 
color charge Tj with the emission of a gluon of color a from parton i. The action of this 
color charge onto the color space is given by 



Cl . . . Cj . . . 



bl...bi...bm) = Sc.bl ■ ■ ■ T^^b, ■ ■ ■ ^Crahm- 



(B.6) 



Explicit expressions for the T^^ and the color charge algebra are given in [45]. The real 
gluon emission in the soft limit is then given by: 



da 



soft 



27r 



-1 dt 



IGvrs^ ' 



(B.7) 



where the phase space integrals Xjj are given in Eq. (|B.2| )^ and Fij denote color correlated 
amplitudes. At O(a^a), Fij are given by 



M 



1,0 

12-s>34 



2 Re 



M 



0,1 ci...bi...bj ...C4 
12->34 



rpa rpa .,l,0ci...Ci...Cj...C4 

-'6iCi-'6jc/^' 12^34 



(B. 



^In this case A is the fictitious gluon mass. 
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In the case of squark-squark production the tree-level amplitude can be decomposed ac- 
cording to their color structure as 



,jl,J C1C2C3C4 
12^34 



Sc,cA2C4M^^ + 5c,cJc2c,M^\ = (1, 0), (0, 1), 

0,1 



Ml'' = \ [M]f - \M}f 



M 



1,0 



1 



M 



u 

1,0 



(B.9) 



where ■M.ijj' are the amplitudes corresponding to the u-channel and t-channel diagrams, 
respectively. In this case the color correlated amplitudes J-ij are given by 



J^i2 — -^34 — 4 

-^13 = -^24 = -12 (A^?'M M{ 



1,0 



A^?'M A^^''^+ (TW^'M M{' 



(1,0 



^0,1 



,1,0 



14 



12 



,1,0 



<o,i 



,1,0 



1,0 



.0,1 



Ml'' 



+ 4 



12 LM 



<o,i 



7W 



1,0 
2 > 



1,0 



^0,1 



(B.IO) 



where in the last case i = 1, . . . , 4. 

The partonic cross section in the collinear region is again given by a convolution inte- 
gral, similar to Eq. (B.4), 



'^'^12^343 



OsCf 



coll 



vr 



1-5,, 



dz Kqq{z,s) d'5'x2_>34 



(zs), 



(B.ll) 



with Ci? = 4/3. zq ^-nd defined in Eq. ( [8.41) . 

B.3 Collinear quark bremsstrahlung 

For real quark emission, only initial-state collinear singularities arise. The partonic cross 
section in the collinear region is given by a convolution integral. 



da 



2,1 



coll 



27r 



'12->34g('^) 

with Tp = 1/2. Kqg is given by 

Kgg = {z^ + (1 - zf) In 
while zq is defined according to Eq 



Az Kqg{z, s)dal2_^^^{zs) 



12->34V^"';i 



(B.12) 



20 



s5g{l - z) 

Iniq 



+ 2z(l-z) 



(B.13) 
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C. Feynman diagrams 



In this Appendix, we list all parton-level Feynman diagrams for the EW contributions to 
the generic process qq' — )• qaq'p with q^q' = {u, d, c, s}. The complete list of LO Feynman 
diagrams is given in Figure |l|. 




Figure 13: Virtual corrections (I): EW one-loop insertions into QCD Born diagrams. Q and Q' 
denote the SU(2) partner of quark q and g', respectively. We use generic labels = J,Z and 
5'° = h^,H^^G^,A^. If the chirality of an internal squark is not specified, it can be any. The 
diagram containing the four-squark vertex has to be taken at O^asO). The diagrams in the third 
line contribute only for u = {u,c}, d = {rf, s}- The last two diagrams contain the counter terms, 
whose rcnormalization constants have to be evaluated at 0{a). For q — q' crossed diagrams have 
to be taken into account. 
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Figure 14: Virtual corrections (II): QCD one-loop insertions into EW Born diagrams. For q~q' 
crossed diagrams have to be taken into account. The diagrams containing the four squark interaction 
vertex have to be evaluated at 0{asOi). The chargino-mcdiated diagrams only contribute for u = 
{u, c}, d = {d, s}. The renormalization constants appearing in the counterterm diagrams (last line) 
have to be evaluated at 0{as)- 
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q Qa q Qa q Qa q qa q g qa 





q'p q q q 



Figure 15: Virtual corrections (III): QCD one-loop insertions into QCD Born diagrams. For 
q = q' crossed diagrams have to be taken into account. Here, Qi can be any of the six quark flavors. 
The diagram containing the four squark vertex has to be evaluated at ©(a^). The renormalization 
constants appearing in the counter term diagrams (last three diagrams) have to be evaluated at 
0{as)-, i.e. the strong sector has to be rcnormalizcd. 




Figure 16: Feynman diagrams for real photon emission. For q ~ q' crossed diagrams have to be 
taken into account. Diagrams in the last row only contribute for u = {u, c}, d = {d, s}. 
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Figure 17: Feynman diagrams for real gluon emission. For q — q' crossed diagrams have to be 
taken into account. Diagrams in the last row only contribute for u = {u, c}, d = {d, s}. 




Figure 18: Feynman diagrams for real quark emission. For q ^ q' diagrams with q and q' exchanged 
have to be considered, too. For q — q' crossed diagrams have to be taken into account. Diagrams 
in the last row only contribute for u = {u, c}, d = {d, s}. 
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